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Abstract. 

Wc discuss a real time, non-perturbative computation of the transverse dynamics 
of gluon fields at central rapidities in very high energy nuclear collisions. 

INTRODUCTION 

This year (1999) the Relativistic Heavy Ion Colhder (RHIC) at Brookhaven will 
begin colliding beams of gold ions at center of mass energies of 200 GeV/nucleon. In 
slightly over 5 years from now, the LHC collider at CERN will do the same at center 
of mass energies of about 5.5 TeV/nucleon. At these energies, the appropriate 
basis to describe the colliding nuclei is that of partons, the quarks and gluons 
that constitute a nucleon, rather than a hadronic basis. An objective of the above 
mentioned experiments is to determine whether the partons, confined in nucleons 
prior to the collision, are liberated after the collision to form fleetingly, in the 
relatively large nuclear volume, an equilibrated state of matter popularly known as 
the quark gluon plasma. 

Clearly, the space-time evolution, and possible equilibration, of matter produced 
in a nuclear collision must depend on the initial conditions. These are given by 
the parton distributions inside each of the nuclei. In perturbative QCD, the "fac- 
torized" expression for the multiplicity or energy distribution of a high pt jet, with 
P)- ^ ^/s, is obtained by convolving the parton distributions in each of the nuclei, 
at the hard scale of interest, with the elementary parton-parton cross section. If 
X = Pt/y/s is not too small, the factorized expression is reliable. With the mea- 
sured nuclear structure functions, one can then compute the multiplicity and energy 
distributions of the jets produced [1]. 
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However, for x « 1 (corresponding to the transverse momentum range 
^QCD « Pt « Vs) the factorization formula for energy and multiphcity dis- 
tributions breaks down. Simply put, this is because partons in one nucleus can 
resolve more than one parton in the other [2]. The parton densities in the nu- 
clei become very large and may even saturate at sufficiently small x. The regime 
of high parton densities is a novel regime in QCD where, although the coupling 
constant may be small, the fields strengths are large enough for the physics to be 
non-perturbative [3]. 

The precise x value at which the above mentioned leading twist factorization 
breaks down is not clear. There are hints from from HERA that parton saturation 
is already seen in the data for x ^ 10~^ and fa 4 GeV^ [4]. If this result is 
robust, similar efi^ects may be seen in nuclei at larger values of x, even x ~ 10"^. 
Their relevance for RHIC (and especially LHC) cannot then be ignored. 

The efi'ccts of high parton densities in the central rapidity region of nuclear 
collisions can be studied in a model which is based on an effective field theory (EFT) 
approach to QCD at small x [5] . The model describes the time evolution of gauge 
fields in a nuclear collision. It takes into account, self-consistently, interference 
effects (which are also responsible for shadowing in deeply inelastic scattering) that 
become important at small x. Another nice feature is that it provides a space-time 
picture of the nuclear collision. This feature would be extremely useful if the gauge 
fields at late times were to provide the initial conditions for a parton cascade [6] or 
for hydrodynamic evolution if it can be determined that the matter produced has 
equilibrated [7]. 

The model is formulated in the infinite momentum frame — > oo and light 
cone gauge A'^ = 0. It contains a dimensionful parameter /x^, defined to be the 
color charge squared per unit area. 



Here q, g stand for the nucleon quark and gluon structure functions at the resolution 
scale Q of the physical process of interest. Also, above xq = Q/^. Using the 
HERA structure function data, Gyulassy and McLerran [8] estimated that < 1 
GeV for LHC energies and < 0.5 GeV at RHIC. Thus a window of applicabihty 
for weak coupling techniques does exist, and higher order calculations will tell us 
if it is smaller or larger than the simple classical estimate. 

An interesting property of the light cone gauge is that final state interactions are 
absent! Kovchegov and Mueller [9] showed that the effects of final state interactions, 
as seen in a covariant gauge computation, are already contained in the nuclear 
wavefunction in fight cone gauge. This non-trivial observation is at the heart of 
the approach described in this talk. Finally, we should alert the reader to alternative 
approaches to the one described here pursued in Refs. [10]- [11]. 
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CLASSICAL MODEL OF GLUON PRODUCTION 



At very high energies, the hard valence quark (and gluon) modes are highly 
Lorentz contracted, static sources of color charge for the wee parton, Weizsacker- 
Williams, modes in the nuclei. The sources are described by the current 

r^^n) = 5''+pl{n)5{x-) + 5^-p^(ri)5(a:+) , (2) 

where pi(2) correspond to the color charge densities of the hard modes in nucleus 
1 (nucleus 2) respectively. The classical field describing the small x modes in the 
EFT is obtained by solving the Yang-Mills equations in the presence of the two 
sources. We have then 

D^F^^ = r . (3) 

The small x glue distribution is simply related to the Fourier transform A^ikt) of 
the solution to the above equation by < A1{kt)A'^^{kt) >p. 

The above averaging over the classical charge distributions is defined by 

(O), = / 0(p, p.) exp (- / a^r.^^MM±^^ . (4) 

We have assumed identical nuclei with equal Gaussian weights g^fi^. 
Before the nuclei collide {t < 0), a solution of the equations of motion is 

A^^O;A'^ 9{x-)9{-x+)a{{rt) + 9{x+)9{-x-)a2{rt) , (5) 

where ag{rt) {q = 1,2 denote the labels of the nuclei and i = 1,2 are the two 
transverse Lorentz indices) are pure gauge fields defined through the gauge trans- 
formation parameters h.q{r],rt) [8] 

Q/*(rt) = -(Pe ■'^•'proj j V I Pe =^''p™J 

Here 1] = itT^proj =F log(a;^/a;|^j.Qj) is the rapidity of the nucleus moving along the 
positive (negative) light cone with the gluon field a\(^2)- '^^^ ^qiv^ ^t) are determined 
by the color charge distributions A^A^ = pq (q=l,2) with being the Laplacian 
in the perpendicular plane. 

For t > the solution is no longer pure gauge. Working in the Schwinger gauge 
= x^A~ + x~A'^ = 0, Kovner, McLerran and Weigert [12] found that with the 
ansatz 

A± = ±x^a{T, rt) ; A' = al(T, n) , (7) 

where r = y/2x'^x~, Eq. 3 could be written more simply in terms of a and a±. Note 
that these fields are independent of 77-the solutions are explicitly boost invariant in 
the forward light cone! 




The initial conditions for the fields a(r, rj) and at r = are obtained by 
matching the equations of motion (Eq. 3) at the point a;^ = and along the 

boundaries x"*" = 0,a;~ > and x~ = 0,a;"'" > 0. Remarkably, for such singular 
sources, there exist a set of non-singular initial conditions for the smooth evolution 
of the classical fields in the forward light cone. One obtains 

i 

a'^_\r=o = a{ + al] a\r=o ^ -[a\, a^] . (8) 

Gyulassy and McLerran have shown [8] that even when the fields a\ 2 before the 
collision are smeared out in rapidity, to properly account for singular contact terms 
in the equations of motion, the above boundary conditions remain unchanged. 
Further, since as mentioned the equations are very singular at r = 0, the only 
condition on the derivatives of the fields that would lead to regular solutions are 

dra\r=0, draWr=0 = 0. 

Perturbative solutions of the Yang-Mills equations to order in the color charge 
density (or equivalently to second order in asfJ'/kt) were found, and at late times, 
after averaging over the Gaussian sources, the number distribution of classical glu- 
ons was found to be [12,8,13,14] 



dyd?kt (27r)4 kt 

where L{kt, A) is an infrared divergent function at the scale A. This result agrees 
with the quantum bremsstrahlung formula of Gunion and Bertsch [15]. Also, Guo 
has shown that the above result is equivalent to the perturbative QCD factorized 
result for the process qq qqg [16]. 

From the above expression, it is clear that distributions are very sensitive to 
L{kt, A). Usually, as in Gunion and Bertsch, this divergence is absorbed in a non- 
perturbative form factor. What is novel about the classical approach is that, at suf- 
ficiently high energies, the non-linearities in the Yang-Mills fields self-consistently 
regulate this infrared divergence. To confirm this claim, one needs to solve the 
Yang-Mills equations to all orders in as/J^/h. A non-perturbative solution of the 
Yang-Mills equations on a two dimensional transverse lattice was performed by 
us [17] and is described below. 



REAL TIME SIMULATIONS OF YANG-MILLS I: 
LATTICE FORMULATION 

We have seen above that the Yang-Mills equations are boost invariant. This is 
a consequence of the sources being 5-functions on the light cone- the nuclei are 
assumed to move with the speed of light! Since this is not the case, boost invariance 
is only approximate. It should, however, be a good assumption at the energies of 
interest-especially at central rapidities. 



The boost invariance assumption simplifies our numerical work considerably. We 
now have a 2+1-dimensional theory and all the dynamics is restricted to the trans- 
verse plane (we assume 1] = 0). The Yang-Mills equations are most conveniently 
solved by fixing = gauge and solving Hamilton's equations. Gauge invariance 
is ensured by defining the theory, in the usual way, on a 2-dimensional transverse 
lattice. The lattice Hamiltonian is the Kogut-Susskind Hamiltonian for gauge fields 
coupled to an adjoint scalar 

Hl^^ E £^r£^r + r E (l - ^Tr[/o) , 

+ i E Tr - t/.,n*,+nt^,,)' + ^ E Trp?, (10) 

j,n 3 

where Ei are generators of right covariant derivatives on the group and Uj^n is a 
component of the usual SU(2) matrices corresponding to a link from the site j in the 
direction n. The first two terms correspond to the contributions to the Hamiltonian 
from the chromoelectric and chromomagnetic field strengths respectively. Also, 
above $ = $"(t" is the adjoint scalar field with its conjugate momentum p = •p°'a°' . 

Lattice equations of motion follow directly from Hj^ of Eq. 10. For any dynamical 
variable v, with no explicit time dependence, v — {Hl, v}, where v is the derivative 
with respect to r, and {} denote Poisson brackets. We take Ei, Ui, pj, and $j 
as independent dynamical variables, whose only nonvanishing Poisson brackets are 
{p,",$5} = 5,j6ab; {E^,Um} = -iSimUia'^; {^^r,^} = 25;„e„beEf (no summing of 
repeated indices) . The equations of motion are consistent with a set of local con- 
straints (Gauss' laws). Their evolution in r after the nuclear collision is determined 
by Hamilton's equations and their values at the initial time r = 0. 

The initial conditions on the lattice are the constraints on the longitudinal gauge 
potential A"^ and the transverse link matrices U± at r = 0. The longitudinal gauge 
potentials are zero outside the light cone and satisfy the Schwinger gauge condition 
^4"^ = inside the light cone x± > 0. They can be written, as in the continuum 
case (see Eq. 7), as 

= ±x^e(x+)e(x-)a(T, xt) . (11) 

The transverse link matrices are, for each nucleus, pure gauges before the collision. 
This fact is reflected by writing 

= e{-x+)e{-x-)i + e{x+)e{x-)u{T) + e{-x+)e{x-)Ui + e{x+)e{-x-)U2 , 

(12) 

where are pure gauge. 

The pure gauges are defined on the lattice as follows. To each lattice site j 
we assign two SU(A^c) matrices Vij and ¥2^. Each of these two defines a pure 
gauge lattice gauge configuration with the link variables = VqjVgj_^_j^ where 



q — 1,2 labels the two nuclei. Also, as in the continuum, the gauge transformation 
matrices V^j are determined by the color charge distribution pgj of the nuclei, 
normally distributed with the standard deviation g^/ij^: 

^[pj°^exp|^-^x:pL)- (13) 

Parametrizing Vgj as exp(iA^) with Hermitean traceless Aj, we then obtain Aj by 
solving the lattice Poisson equation 

= E {^U + Aj-n - 2Aj) = p,j. (14) 

n 

The correct continuum solution (Eqs. 5 and 7) for the transverse fields A± is re- 
covered by taking the formal continuum limit of Eq. 12. 

Using the general representation of the gauge fields in Eqs. 11 and 12, we shall 
now state the initial conditions for them at r = 0. 

U±\r=o = (C/i + U2)iUl + Ul)-' ; Eii^o = 0. 
Pj\r=o = 2q; ; $j = , (15) 

where U± is defined as exp{ia±a±). The above initial conditions are obtained by 
matching the lattice equations of motion in the four light cone regions at r = 0. 
For details we refer the reader to Ref. [17]. 



REAL TIME SIMULATIONS OF YANG-MILLS II: 
DISCUSSION OF RESULTS 

In Ref. [17], we reported results of simulations of the time evolution of classical 
fields in a 2-|-l-dimensional SU(2) gauge theory described by the Hamiltonian in 
Eq. 10. The simulations were carried out on transverse lattices ranging from 20 x 20 
sites to 160 X 160 sites. The lattice results depend on one dimensionlcss parameter, 
g^nL [18]. The parameter L corresponds to the size of the nucleus; ii defined in 
Eq. (1) is determined by the size of the nucleus, the energy of the nucleus, and the 
hard scale Q of interest; the strong coupling g runs as a function of either Q or /i 
depending on which is greater. 

Thus, once given the energy and size of the incoming nuclei, and the hard scale 
of interest, we can determine the evolution of gauge fields in the central rapidity 
region. Consider, for instance, coUiding two gold nuclei (A~ 200) at RHIC and 
LHC energies. Approximating = vri?^, we obtain L = 11.6 fm. For the hard 
scales of interest, Q ^ 1-2 GeV, /j, - 0.5 GeV at RHIC (/^ - 1 GeV at LHC). One 
then obtains 



g'^pL ^ 120 (RHIC) 
fti240 (LHC). 



(16) 



Above we have chosen = 2 (or equivalently, as ~ 0.3). 

On the lattice, the lattice couphng is g'^{iia){L/ a) = q'^/jllN. The continuum 
limit is obtained by keeping g'^fi^N fixed (to the physical value of interest-as in 
Eq. 16) and taking fiL zero. It appears from our simulations that we are in the 
weak coupling regime for /il = 0.017, 0.035 in lattice units. For the physical values 
of g^/iL above, these would correspond to lattices an order of magnitude larger 
than those considered so far. Detailed simulations on the above physical scenario 
will be reported at a later date. 

We now turn to an issue of some concern; whether quantities of interest have a 
continuum limit (in the above defined sense) in the classical theory. For instance, 
in thermal field theories, it is not clear that dynamical quantities such as auto- 
correlation functions have a well defined limit as the lattice spacing a — > 0. In the 
EFT described here, there is reason to be more optimistic. 

Consider the following gauge invariant quantity; the energy density — 
Ek/N'^ of the scalar field on the lattice (in units of nj). It is plotted as a function 
of the lattice size N (in units of the lattice spacing) in Fig. 1 for hl — 0.0177, 0.035. 




FIGURE 1. The lattice size dependence of the scalar kinetic energy density, expressed in units 
of /i^ for fj, = 0.0177 (pluses) and fj, = 0.035 (diamonds). The solid line is the LPTh prediction. 
The error bars are smaller than the plotting symbols. 



The solid line is the prediction from lattice perturbation theory (LPTh.). It is given 

by 



pV = 6 ) E 



E 



sin(/^) sin(/^/ ) 
A^(0 



+ 16 E 



sin2(^)sin2(^ 



A2(0 



, (17) 



where In — 27rn/N and A(Z) = 2X^„=i 2(l — cos(Z„)) is the usual lattice Laplacian. 
The continuum limit of the above equation has the form p'^p"' — > A + B\og\L/a), 
where A and B are constants that can be determined from the above equation. 



Prom Fig. 1, it can be seen that the kinetic energy density diverges as log^(L/a)- 
i.e., it does not have a well-defined continuum limit. This divergence is softer 

than in a thermal theory, where the energy density diverges as an inverse power of 
the lattice spacing. If wc subtract the perturbative result in Eq. 17 at some scale 
^nanperu the resulting expression will have a continuum limit but will depend on 
^ncmpert- The trick is to find a h.ncm,pert for which the contribution to the "non- 
perturbative" scalar kinetic energy density converges as /i^ — > 0. In Table. 1, 
we show the results, at r = 0, of simulations where g^^L = 33.6 is kept fixed 
(note: here we fix 51 = 1). The lattice cut-offs (recall that kt = 27Tn/L) of ricut — 
10,15,20 correspond to A.nonpert = 1.08,1.62,2.16 GeV respectively. For ricut — 

TABLE 1. Scalar k.e. density ]v%^ A*^ 
for fixed g^fiL = 33.6. Columns 3-5 are the 
scalar. ke.d after subtracting LPTh. contribu- 
tion for n > ricut in Eq. 17. 



Ml 




ncut=W 


ricut =15 


ncut='20 


.84 


07.24C-02 


.0341 


.0562 


.0702 


.56 


10.20e-02 


.0407 


.0636 


.0791 


.42 


12.38e-02 


.0426 


.0694 


.0852 


.21 


17.60e-02 


.0268 


.0702 


.0945 



15, the non-perturbative contribution to the scalar. ke. density appears to converge 
to a constant value as hl is decreased, keeping g^/iL fixed. This result must of 
course be confirmed for larger lattices and other values of g'^fiL. Caveat aside, our 
results appear to suggest that for a particular Anonpert there is a non-perturbative 
contribution to the energy density that survives in the continuum limit. 

Presumably, the subtraction described above may also be performed for static 
quantities in a thermal field theory. It is unlikely that this procedure there is 
reliable for dynamic quantities. In our case, the scalar kinetic energy is a dynamic 
quantity evolving in time. If the above mentioned procedure is to be use, it must 
also be valid at later times. We are optimistic that this may be the case because 
our simulations suggest that the time evolution of hard modes decouples from that 
of the soft modes. This is illustrated very clearly by Fig. 2 of the most recent of 
our papers cited in Ref. [17]. 

The parameter Anonpert niay be given a physical interpretation. At sufficiently 
high energies, one may be able to relate it to the physical scale at which one begins 
to see deviations from perturbative predictions (note: at very high energies one 
expects Anonpert >> Aqcd) duc to high parton density effects. 

We are currently studying whether the non-perturbative contribution to the 
scalar energy density also has a robust continuum limit at late times [19]. 
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